Diffusion-limited reactions in crowded environments 
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Diffusion-limited reactions (DLR) are usually described within the Smoluchowski theory, which 
neglects interactions between the diffusing components. We propose a first extension of such frame- 
work that incorporates excluded-volume effects, considering hard diffusing agents in the presence 
of an absorbing sink. For large values of the absorber-to-diffuser size ratio Rs, the encounter rate 
obtained from the simulations is in very good agreement with a simple generalization of the Smolu- 
chowski equation at high densities. Reducing Rs, the rate gets substantially depressed at interme- 
diate packing and become even non-monotonic for Rs <^ 1. Concurrently with the saturation of the 
rate, stationary density waves set in close to the absorbing center. A mean-field, effective-potential 
analysis of the Smoluchowski equation under crowding sheds light on the subtle link between such 
ordering and the slowing down of the encounter dynamics. Finally, we show how an infinitesimal 
amount of non-reacting impurities can equally slow down dramatically the reaction. 
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Processes such as gelation, coagulation, crystallization 
or self-assembly in colloidal or polymer systems, thin-film 
growth in materials science and chemical reactions in bi- 
ology involve, as a first step, the random movement of 
particles suspended in a fluid [1]. The components dif- 
fuse and when they come into contact, a reaction may be 
triggered, such as aggregation or assembling. As an ex- 
ample, the usual paradigm for biochemical reactions as- 
sumes the formation of an encounter complex, which may 
subsequently undergo chemical transformation, yielding 
the product. In most cases, the fixation step proceeds 
much faster than the encounter, in which case the reac- 
tion is said to be diffusion-limited. Such reactions are 
commonly found in biochemical processes such as en- 
zyme catalysis, protein aggregation and complexation in 
cells [2]. 

The simplest model of diffusion-limited encounter in 
three dimensions has been formulated almost one cen- 
tury ago by Smoluchowski under the hypothesis of non- 
interacting, spherical and chemically isotropic reactants, 
a diffusing particle (p) and a sink (s) [3|, leading to the 
Smoluchowski absorption rate 
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where pao is the (relative) bulk density of the reactants, 
D = Ds + Dp is the coefficient of relative diffusion and 
R = Rs + Rp the encounter distance. 

The Smoluchowski theory is still nowadays the main 
theoretical framework within which the aforementioned 
processes are analyzed. However, this approach is strictly 
valid only for ideal, infinitely diluted solutions, while 
most systems become interesting at concentrations far 
from the ideal gas limit. Cells, for instance, contain a 
large number of proteins, nucleic acids, and other smaller 



molecules that occupy up to 30-40% of the available vol- 
ume |1] and that cannot overlap with each other. As a 
matter of fact, crowding effects are expected to impact 
profoundly on the thermodynamics and kinetics of bi- 
ological processes in vivo [IHS], such as protein folding 
and stability [7] and aggregation [5]. Yet, the effects of 
crowding on diffusion-limited processes have been exam- 
ined only in the case of low density of diffusing parti- 
cles [3] or for tracers diffusing in a medium of inert parti- 
cles [T0Hl3] . In colloidal science, instead, the properties 
of crowded environments have been investigated widely. 
Experiments on concentrated suspensions of nearly hard- 
core particles revealed unexpected phenomena both in 
their thermodynamics and in the dynamics, spotlighting 
a subtle entwining between packing, structuration and 
dynamics |14j . 

In this paper, we generalize the classic Smoluchowski 
problem to arbitrary crowding conditions by means of a 
novel computational scheme that allows to efficiently ex- 
plore the effects of increasing packing on the encounter 
dynamics. Typical numerical simulations of diffusion- 
limited reactions employ different declensions of the 
Brownian dynamics (BD) algorithm [121 HH]- Here, we 
consider a liquid composed of hard spheres of radius 
Rp described through Event Driven Brownian Dynam- 
ics {EDBD) [T7HT9] . This technique allows to efficiently 
simulate the stochastic dynamics of many-body hard-core 
systems. In order to simulate encounter reactions at fi- 
nite densities, we adapted the EDBD to the configuration 
of an absorbing sink of radius Rs located at the center 
of a spherical bounding box with constant-flux boundary 
conditions (Fig. [T]) [20^ . In this way, by tuning the density 
in the box, we can measure the encounter rate n and the 
stationary density profile p{r), from the infinite-dilution 
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FIG. 1: (Color online) Simulation scheme. The reacting par- 
ticles (yellow spheres) diffuse in a spherical box. Upon reach- 
ing the absorbing sink (red sphere) located at the center of 
the box, they get absorbed and are reinserted in the external 
buffer (green sphere) . 
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Smoluchowski limit, namely 



/^S, PS{A = Poo [1 - {Rs + Rp)IA 



(2) 



up to high crowding conditions. In the following, densi- 
ties are expressed in terms of the corresponding packing 
fraction of the spheres 0, and lengths in units of Kp. 

As a first investigation of crowding effects, we sim- 
ulated the case of a large sink-to-particle diameter ratio 
(i?s = 8/3). At low densities (0 < 0.1), the rate increases 
linearly, in agreement with the BD and analytical results 
by Dzubiella and McCammon [5], following the first or- 
der virial correction to the Smoluchowski rate. Increas- 
ing the crowding further, the encounter rate is strongly 
enhanced, reaching even an eight-fold magnification at 
(/) = 0.4 with respect to the infinitely diluted limit. 

In order to account for the strongly nonlinear increase 
of the rate at intermediate crowding, one can consider 
the Smoluchowski equation and take inter-particle inter- 
actions into account in the expression of the net flux J 
through a generalized Stokes-Einstein relation for the col- 
lective diffusion coefficient [21]. At first order, crowd- 
ing effects can be accounted for through the density de- 
pendence of the pressure in the bulk liquid n(p) (van- 
ishing at zero density), yielding a generalized diffusiv- 
ity D{p) = i:>orf[/3n(p)]/d/9, where ^'^ = UbT and 
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D{p). The steady-state Smoluchowski 



equation in spherical coordinates reads then 
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with the boundary conditions p{Rs + Rp) 
and limj._>oo p(r) = p^o- Integrating the equation 
that defines the encounter rate, namely K/(47rr^) — 
J{r), in the interval [i?, +oo) with J{r) = \ J\ ~ 
l3Do(dp{r)/dr){dIl{p)/dp) one obtains 
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FIG. 2: (Color online) Encounter rate versus crowding. For 
large sink-to- particle diameter ratios {Rs = 8/3), Eq. Q ac- 
counts for the EDBD results. Reducing Rs, the rate saturates 
at high packing, while it displays a maximum for Rs <^ 1. The 
virial solution and the BD results of Ref. i9j are also shown. 
In the inset, we plot the rates versus Rs for (j) = 0.33 and 
(j) = 0.44, together with the corresponding theoretical limits 
predicted by Eq. Q. 



where I31l{poo)/poo — Z{(f>) is the compressibility fac- 
tor. The latter can be well approximated through the 
Carnahan-Starling (CS) equation of state Z{(f)) = (1 -|- 
4> + (1)'^ - (f>^)/{l - (j>)^ [22^. At moderate densities, the 
virial expansion /3Il{p) — p(l + B2P + 0{p^)) yields the 
Smoluchowski rate at zeroth order and the first-order cor- 
rection k/ Ks = 1 + B2P00 derived in Ref. [SI 

For large values of the sink size, the rate increases 
monotonically with the density in good agreement with 
Eq. Q, which provides an accurate prediction in the 
limit Rs ^ Qo. This is clearly illustrated by inspect- 
ing how the rate varies with R^ at fixed packing fraction 
(inset in Fig. [2|. 

The density gradient between the sink and the bulk 
regions induces a net force towards the sink, the lat- 
ter behaving as a depression in an high pressure envi- 
ronment. The rate is thus directly related to the dif- 
ference in pressure between the diluted region near the 
sink and the denser bulk. Remarkably, however, when 
Rs is decreased the rate gets substantially depressed at 
intermediate crowding and becomes even non-monotonic 
for i?s <C 1, featuring a maximum and a drop at high 
crowding. 

The stationary density distributions of the diffusing 
particles around the absorbing boundary surface offer a 
key to understand the drop in the rate observed at in- 
termediate packing when reducing Rg. As illustrated in 
Fig. [3| we found that (i) oscillations set in progressively 
while increasing 4> and (ii) the structuration is strongly 
enhanced for small values of Rs- More specifically, the 
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rate reads 



FIG. 3: (Color online) Structuration of the density profiles 
under crowding. The profiles are plotted at increasing bulk 
packing fraction cj> for Rs = 2/3 (left) and Rs — 1/30 (right). 
Structuration is enhanced in the case of the smaller sink, but 
the liquid gets structured at approximately the same distance 
A from the contact distance Rs + Rp in the two cases. The 
density oscillations are in phase with the radial distribution 
function g{r) of a hard-sphere liquid at intermediate densities. 



accumulation of particles around the sink corresponding 
to the first maximum in the density profiles occurs at a 
distance A' ~ {Rg + Rp) + A, where A does not depend 
on Rg. 

As it is clearly shown in Fig. l3J the progressive struc- 
turation of the system bears the clear signature of the 
radial distribution function g{r) of a hard-sphere liquid 
at intermediate densities. These effect is by no mean triv- 
ial, as g(r) provides a measure of spatial correlations in 
the reference frame of a moving particle. Here, transla- 
tional symmetry is broken by the sink, which causes the 
emergence of denser layer near the absorber, triggering 
the onset of stationary density waves. 

The simultaneous onset of density oscillations and rate 
falloff observed upon reducing Rg suggests that the or- 
ganized accumulation of particles is the consequence of 
a competition among the particles diffusing in the prox- 
imity of the sink. To this regard, it is important to note 
that, A being a constant, the particles compete for the 
sink at a distance from the absorber that depends on the 
size of the diffusers only, and not on the size of the latter. 

In order to emphasize the connection between the in- 
creasing structuration of the liquid and the slowing down 
of the encounter process it is possible to regard the onset 
of density fluctuations as marking the appearance of an 
effective potential between the sink and individual par- 
ticles. Within such representation, the accumulation of 
particles at A will act on a freely diffusing particle as 
an energy barrier V{r) — f{r;(/),Rs) (Fig. kl(a)), whose 
height should depend parametrically on both the bulk 
concentration and the sink size Rg to account for the 
increasing structuration of the density profile for large 
values of and small values of Rg- In general, in the pres- 
ence of an interaction potential V{r), the Smoluchowski 
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It is possible to generalize Eq. (Is]) in order to account 
for the effects of crowding by noting that finite densities 
imply K ex /3n(/9oo)j as shown by Eq. ([4]), instead of k oc 
Poo, as from Eq. ([I]) (free particles) and Eq. ([5| (sink- 
particle potential). As a consequence, we may write 
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. It is reasonable to assume that the height V of the 
potential energy barrier reflects the excess pressure that 
accumulates in the vicinity of the sink with respect to 
the ideal case as the packing fraction increases. As a 
consequence, we assume 
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where rj is related to the local density at a distance A 
(Fig. |4] (b)) and Ilid{p) = Pp is the ideal gas pressure. 
For simplicity, we take rj = P5(A) which introduces the 
dependence on Rg in the effective potential. 

Remarkably, a one-parameter fit of Eq. (l6| to the 
numerical data accounts extremely well for both the (p 
and Rg dependence of the rate (Fig. l4] c) with a sin- 
gle best-fit value of the floating parameter (c = 2.1) for 
< < 0.52 and 1/30 < -R^ < 8/3, thus substantiat- 
ing our theoretical analysis. Despite its simplicity, this 
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FIG. 4: (Color online) Effective potential description of DLRs 
under crowding. The effective potential V (</!>, Rs) (a) accounts 
for the accumulation of particles at distance A, resulting in 
the oscillations in the density profiles (b). The solid lines 
in panel (b) are plots of the Smoluchowski profiles ([T|, em- 
phasizing the Rs dependence of the density at A. (c) The 
one-parameter fits obtained through Eq. Q are in excellent 
agreement with the numeric results for all values of -Rs, pro- 
viding a quantitative link between the drop in the rate and 
the onset of stationary density oscillations. 
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open issue. Our novel approach clearly highlights the 
physical ingredients necessary for the derivation of a 
self-consistent Smolucho-wski theory valid for cro-wded 
environments. Moreover, -while jamming phenomena like 
arches at bottleneck or at escape issues are -well kno-wn 
in granular media or traffic studies [57], this is, to the 
best of our knowledge, the first time that a complex 
ordering in the presence of an absorbing boundary is 
reported for a system that is still in a liquid-like state. 
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FIG. 5: (Color online) Encounter rate versus fraction of non- 
reacting impurities x for Rs — 2/3. Inset: the density profile 
for impurities (;/>/(r) for x = 1/2000. 



approach seems to capture the essence of the problem, 
namely a density-weighted competition among the parti- 
cles diffusing around the sink. Diffusion-limited reactions 
in crowded environments are thus the result of a subtle 
combination of the liquid properties of the diffusers with 
the underlying long-range Smoluchowski description. 

As a first extension of our numerical scheme to more 
complicated diffusion- limited reactions [U [HI [24], we 
investigated the effects of introducing a finite amount of 
non-reacting hard sphere impurities (/) in the system. 
We found that an infinitesimal fraction {x = 1/2000) is 
sufficient to slow down considerably the reaction with 
respect to the one-component case (Fig. [s]). While 
at low crowding the rate and the density profiles are 
insensitive to the presence of non-reacting particles, an 
accumulation of the latter against the sink takes place 
when increasing above 0.35. As a result, the sink 
gets screened to the reacting particles and the rate falls 
off. Indeed, when (j) is increased, the impurities get 
trapped at the sink surface because of the large density 
gradient, and need to overcome a substantial energy 
barrier in order to get back to the bulk region |5S]. The 
extreme sensitivity of the reaction rate to infinitesimal 
amounts of impurities might for instance be of great 
interest in the context of molecular trafficking across 
the nuclear pore complex, a protein complex responsible 
for protected exchange between the nucleus and the 
cytoplasm, or for preventing the transport of undesirable 
material through the nuclear envelope [SB] . 

In conclusion, in this letter we tackled a fundamental 



[1] F. Family, D.P. Landau Kinetics of Aggregation and 

Gelation ( North- Holland, The Netherlands, 1984). 
[2] S.A. Rice, Diffusion-Limited Reactions. (Elsevier, 1985). 
[3] M. V. Smoluchowski, Z. Rhys. Chem. 92, 129 (1916). 
[4] R.J. Ellis, A.P. Minton, Nature 425, 27 (2003). 
[5] T.A. Ryan et al. Science 239, 61 (1988). 
[6] S.B. Zimmerman, A.P. Minton, Annu Rev Biophys 

Biomol Struct 22, 27 (1993). 
[7] M. S. Cheung, D. Khmov, D. Thirumalai, Rroc Natl Acad 

Sci 102, 4753 (2005). 
[8] A. P. Minton, Curr Opm Struct Biol, 10, 34 (2000). 
[9] J. Dzubiella, J. A. McCammon, J. Chem. Rhys. 122, 
184902 (2005). 
[10] W. Dong, F. Baros, J. C. Andre, J. Chem. Rhys. 91, 

4643 (1989). 
[11] J. Sun, H. Weinstein, J. Chem. Rhys. 127, 155105 (2007). 
[12] J.D. Schmit, E. Kamber, J. Kondev, Rhys. Rev. Lett. 

102, 218302 (2009). 
[13] J.S. Kim, A. Yethiraj, Biophys. J. 98, 951 (2010). 
[14] P.N. Pusey, W. van Megen, Nature 320, 340 (1986). 
[15] D.L. Ermak, J.A. McCammon, J. Chem. Rhys. 69, 1352 

(1978). 
[16] B. Cichocki, K. Hinsen, Rhys. A 187, 133 (1992). 
[17] G. Foffi et al., Rhys. Rev. Lett. 94, 078301 (2005). 
[18] A. Scala, T. Voigtmann, C. De Michele, J. Chem. Rhys. 

126, 134109 (2007). 
[19] C. De Michele J. Comput. Rhys. 229, 3276 (2010). 
[20] N. Dorsaz et al. J. Rhys. Cond. Mat. 22, 104116 (2010). 
[21] J. Dhont, An Introduction to Dynamics of Colloids. (El- 
sevier, Amsterdam, 1996). 
[22] N. Carnahan, K. Starling, J. Chem. Rhys. 51, 635 (1969). 
[23] P. L. Houston. Chemical Kinetics and Reaction Dynam- 
ics. Dover Publications, 1985. 
[24] A. Zaccone et al. Rhys. Rev. E 80, 051404 (2010). 
[25] Details of the simulations reported in a forthcoming 

study. 
[26] F. Alber et al.. Nature 450, 695 (2007) 
[27] D. Helbing, I. Farkas, T. Vicsek, Nature 407, 487 (2000). 



